Interaction model for the gap equation 
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We explain a form for the rainbow-ladder kernel whose momentum-dependence is consonant with 
modern DSE- and lattice-QCD results, and assess its capability as a tool in hadron physics. In 
every respect tested, this form produces results for observables that are at least equal to the best 
otherwise obtained in a comparable approach. Moreover, it enables the natural extraction of a 
monotonic running-coupling and -gluon-mass. 
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Solving quantum chromodynamics (QCD) presents a 
fundamental problem: never before have we been con- 
fronted by a theory whose elementary excitations are 
not those degrees-of-frcedom readily accessible via ex- 
periment; i.e., are confined. Moreover, there are reasons 
to believe that QCD generates forces which are so strong 
that < 2% of a nucleon's mass can be attributed to the 
current-quark masses that appear in QCD's Lagrangian; 
viz., forces capable of generating mass from nothing, a 
phenomenon known as dynamical chiral symmetry break- 
ing (DCSB). Neither confinement nor DCSB is apparent 
in QCD's Lagrangian. Yet they are dominant in deter- 
mining the observable characteristics of real- world QCD. 
The physics of hadrons is ruled by emergent phenomena, 
such as these, which can only be explained through the 
methods of nonperturbative quantum field theory. 

Confinement and DCSB arc long-distance phenomena. 
Their understanding requires elucidation of the infrared 
behaviour of the interaction between quarks and gluons. 
Much is misapprehended about confinement. It is there- 
fore important to state that the static potential measured 
in simulations of lattice-QCD is not related in any known 
way to the question of light-quark confinement. Light- 
quark creation and annihilation effects arc fundamentally 
nonperturbative in QCD. Hence it is impossible in prin- 
ciple to compute a potential between two light quarks. 
On the other hand, confinement can be related to the 
analytic properties of QCD's Schwinger functions [H-dj], 
so the question of light-quark confinement may be trans- 
lated into the challenge of charting the infrared behavior 
of QCD's ^-function. 

The gap equation is a starting point for analyses of 
DCSB and confinement @]; and it is the basis for formu- 
lating a symmetry-preserving Poincare-covariant treat- 
ment of bound states through Bethe-Salpeter and Fad- 
deev equations [H, [BJ . The gap equation can be written 



Sip)- 1 = Z 2 (i 1 -p + m hm ) 



where: D^ v is the gluon propagator; F„, the quark-gluon 
vertex; , a symbol that represents a Poincare invariant 
regularization of the four-dimensional integral, with A 
the regularization mass-scale; m bm (A), the current-quark 
bare mass; and Zi^iC 2 , A 2 ), respectively, the vertex and 
quark wave function renormalisation constants, with £ 
the renormalisation point. The gap equation's solution 
is the dressed-quark propagator, 



S(p) 
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wherein the mass function, M(p 2 ), is independent of the 
renormalisation point. (We use a Euclidean metric jf|.) 

In realistic studies the model input is expressed in a 
statement about the nature of the gap equation's kernel 
at infrared momenta, since the behaviour at momenta 
k 2 > 2 GcV 2 is fixed by perturbation theory and the 
renormalisation group [7|, |8[. In rainbow-ladder (RL) 
truncation, which is leading-order in the most widely 
used scheme [91. llOj. this amounts to writing (fc = p — q) 



Z ig 2 D^(k)T u (q,p) = k 2 gik 2 )D^ c ik) ly 
= [k 2 Gmik 2 ) + 4ira pClCD ik 2 )] D^ik) lv 



(3) 



9 2 D^ip - q) —1v,S{q) —T v {q,p), 



(1) 



wherein £)J^ e (fc) is the Landau-gauge free-gauge-boson 
propagator; a p QCD(fc 2 ) is a bounded, monotonically- 
decreasing regular continuation of the perturbative-QCD 
running coupling to all values of spacelike- k 2 \ and 
Gmik 2 ) is an Ansatz for the interaction at infrared mo- 
menta: Gm{k 2 ) < 5 pQ cD(fc 2 ) Vfc 2 > 2 GcV 2 . The 
form of Gmik 2 ) determines whether confinement and/or 
DCSB are realised in solutions of the gap equation. (Lan- 
dau gauge is used for many reasons for example, it 
is: a fixed point of the renormalisation group; that gauge 
for which sensitivity to model-dependent differences be- 
tween Ansdtze for the fermion-gauge-boson vertex are 
least noticeable; and a covariant gauge, which is read- 
ily implemented in numerical simulations of lattice regu- 
larised QCD [H[.) 
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The capacity of the Dyson-Schwinger equations 
(DSEs) to unify the explanation of a wide range of me- 
son and baryon observables entails that the pointwise 
behaviour of Gm.(k 2 ) can be constrained through feed- 
back between experiment and theory. This is a practical 
means by which to develop insight into the momentum- 
dependence of QCD's /3-function [Til . 

Following work on confinement |14[ Il5j . the interac- 
tion at small-fc 2 has often been expressed as either an 
integrable infrared singularity, typically Gm(h?) °c S 4 (k), 
or a finite-width approximation to it [jj, la. Il6l - ll8j . The 
following approximation was used in Ref . [7J : 



(4) 



whereas the 5 4 -function itself was used in Ref. [l6[. Nei- 
ther study directly sampled the solutions of the DSEs 
at complex values of their arguments. This changed with 
Ref. Q , which therefore required better control of numer- 
ical procedures and hence employed an equally weighted 
combination of the (5 4 -function and the following finite 
width representation: 



v ' 2tt 2 lo 6 
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The material difference between this form and Eq. ((4]) is 
the inclusion of a multiplicative factor of k 2 . It was in- 
troduced solely in order to tame singularities encountered 
in the numerical treatment of the transverse projection 
operator k 2 D f ™ e (k). 

Desiring additional simplifications in numerical analy- 
sis, the (5 4 -function component of Gm{k 2 ) was completely 
eliminated in Ref. [HI , leaving the infrared behaviour to 
be described by Eq. ([5]) alone; viz., with s = k 2 , 
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T(s), (6) 



where: j m = 12/(33-2JV», N f = 4, A QC d = 0.234 GeV; 
t = c 2 — 1; and J-(s) = {1 — cxp(— s/[4m 2 ])}/s, m t = 
0.5 GeV. This form of the interaction preserves the one- 
loop renormalisation-group behavior of QCD in the gap 
equation; and has since been employed extensively in the 
successful prediction and explanation of hadron observ- 
ables; e.g., Refs. HElll. 

There is an aspect of the interaction in Eq. ([6]) that 
is usually ignored; namely, it produces a kernel for the 
gap equation which possesses a zero at a small timclikc 
value of k 2 , and rapidly becomes very large and negative 
as the magnitude of the timclikc momentum is increased. 
For example, with typical values of the model parameters 
(Duj = (0.72 GeV) 3 , w = 0.4GeV)E 



000 



Eq. 



0fors = -(0.046 GeV) 2 



(7) 



the magnitude of G(s < 0) exceeds its largest spacelike 
value at s = -(0.22 GeV) 2 ; and \Q{s < 0) < 0| grows 
faster than exponentially with decreasing s. 

These facets of the behaviour produced by Eq. (J5|) are 
in stark conflict with the results of modern DSE and 
lattice studies; viz., the gluon propagator is a bounded, 
regular function of spacelike momenta, which achieves 
its maximum value on this domain at k 2 = (30I4321 ] . 
and the dressed-quark-gluon vertex does not possess any 
structure which can qualitatively alter this behaviour [33l . 
l34j . It is thus long overdue to reconsider a functional 
form whose sole raison d'etre was numerical expediency. 

We therefore choose to explore the capacity of 
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as a tool to compute and connect hadron observables. 
This is readily done now owing to improved numerical 
methods for coping with DSE solutions at complex values 
of their arguments [35|. NB. Equation (|8]) cannot be 



expressed via a non- negative spectral density Q . 

Regarding renormalisation of the gap equation, we fol- 
low precisely the procedures of Refs. @, [l8( and use the 
same renormalisation point; i.e., £ = 19 GeV. With gap 
equation solutions in hand for various quark flavours, one 
can solve homogeneous Bethe-Salpeter equations (BSEs) 
for meson Bcthc-Salpctcr amplitudes and therefrom com- 
pute observable meson properties. For example, in the 
isospin-symmctric limit the pion Bcthc-Salpctcr ampli- 
tude is obtained vi£0 



T w (k;P) 



G((k- q y)(k- q yD^(k- q ) 
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where S(£) is the u = <i-quark propagator, P 2 = — m 2 , k 
is the relative momentum between the constituents, and 
one can choose q± = q ± P/2 without loss of generality 
in a Poincare covariant approach. This form of the BSE 
is symmetry-consistent with the gap equation obtained 
through Eq. © [13, H3 ■ All Bethe-Salpeter amplitudes 
are normalised canonically (see, e.g., Eq. (27) in Ref. @j). 

In Table |T] we list computed results for ground-state 
pscudoscalar- and vector-mesons. The meson masses are 
obtained in solving the BSEs. Valid formulae for the 
other quantities, all of which depend linearly on the me- 
son Bethe-Salpeter amplitudes, are presented in Refs. [1, 
I391 . |40| . (NB. The products f^p^ and JkPk describe 
in-pion and in-kaon condensates [1, El|.) The results 
show that observable properties of vector- and flavoured- 
pscudoscalar-mesons computed with Eq. © are practi- 
cally insensitive to variations of 10 G [0.4, 0.6] GeV so long 



Predictions for numerous pseudoscalar- and vector-meson ob- 
servables are approximately independent of uj on the domain 
u) £ [0.3,0.5] so long as one maintains Dui = constant [29ll , 



2 We actually include a factor of 1/ Z| on the left-hand-side of both 
Eqs. lO and JSJl, which additional improvement ensures multi- 
plicative rcnormalisability in solutions of the gap and Bethe- 
Salpctcr equations |36H . 
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(Duj) 1 ' 3 


0.72 


0.8 


0.8 


0.8 


0.8 


UJ 


0.4 


0.4 


0.5 


0.6 


0.7 


m i,d 


0.0037 


0.0034 


0.0034 


0.0034 


0.0034 


mi 


0.084 


0.082 


0.082 


0.082 


0.082 


A(0) 


1.58 


2.07 


1.70 


1.38 


1.16 


M(0) 


0.50 


0.62 


0.52 


0.42 


0.29 




0.138 


0.139 


0.134 


0.136 


0.139 


/« 


0.093 


0.094 


0.093 


0.090 


0.081 


1/2 


0.48 


0.49 


0.49 


0.49 


0.48 




0.496 


0.496 


0.495 


0.497 


0.503 




0.11 


0.11 


0.11 


0.11 


0.10 


1/2 

Pk 


0.54 


0.55 


0.55 


0.55 


0.55 


m p 


0.74 


0.76 


0.74 


0.72 


0.67 


f P 


0.15 


0.14 


0.15 


0.14 


0.12 




1.07 


1.09 


1.08 


1.07 


1.05 


u 


0.18 


0.19 


0.19 


0.19 


0.18 




0.67 


0.67 


0.65 


0.59 


0.46 


1/2 
Per' 


0.52 


0.53 


0.53 


0.51 


0.48 



TABLE I. Results obtained using Eq. ©, compared with one 
representative set computed using Eq. ©. The current-quark 
masses at the renormalisation point C, = 0.19 GeV were fixed 
by requiring a good description of m-K,K- Dimensioned quan- 
tities are reported in GeV. NB. The "<r" listed here is not 
directly comparable with the lightest scalar in the hadron 
spectrum because the rainbow-ladder truncation is a priori 
known to be a poor approximation in this channel [37l . l38l ] . 



as Duo = constant. Furthermore, that there is no reason 
to prefer Eq. (j^J over Eq. ((8|) . 

However, there is reason to prefer Eq. (jHJ over Eq. ([6]). 
Namely, its pointwisc behaviour accords qualitatively 
with results of modern DSE and lattice studies, and it 
can readily be parametrised as follows [3l|, HH 



s(fc 2 ) 



47TQ; flL (fc 2 ) 
k 2 + m 2 Jk 2 ) ' 



ml{k 2 ) 



M 4 
g 



(10) 



with the functions obtained in this way illustrated in 
Fig.[U As one would expect, the infrared scale for 
the running gluon mass increases with increasing u>: 

0.81 GeV for 



= 0.5 GeV; M g 
0.6 GeV. The values of M g are typical [H EH 



M g = 0.67 GeV for uj 

UJ 

Equally naturally, the infrared value of the cou- 
pling is a decreasing function of uj: a R j,{Q)/-K = 15, 
a RL (M 2 )/ir = 3.8 for ui = 0.5 GeV; and a RL {0)/Tr = 9, 
a RL (M 2 )/TT = 2.2 for uj = 0.6 GeV. A context for the 
infrared value of the running coupling required to de- 
scribe meson obscrvablcs in rainbow-ladder truncation 
is readily provided. With nonperturbatively-massless 
gauge bosons, the coupling below which DCSB breaking 
is impossible via the ga p equations in QED and QCD is 
a c /n « 1/3 [361. l42[ I43I] . In a symmetry-preserving reg- 
ularisation of a vector- vector contact-interaction used in 
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FIG. 1. Upper panel - Rainbow-ladder gluon running-mass 
inferred from Eq. © via Eq. (|10[1 . Lower panel - rainbow- 
ladder effective running-coupling inferred from Eq. ©. In 
both panels: ui = 0.5 GeV (solid curve); uj — 0.6 GeV (dashed 
curve) . 



rainbow-ladder truncation, a c /n 0.4 and a description 
of hadron phenomena requires a/ir 1 [44|. With non- 
perturbatively massive gluons and quarks, whose masses 
and couplings run, the infrared strength required to de- 
scribe hadron phenomena in rainbow-ladder truncation 
is unsurprisingly a little larger. Moreover, whilst a di- 
rect comparison between a R r, and a coupling, aQLai, in- 
ferred from quenched-lattice results is not sensible, it is 
nonetheless interesting that aQLat(0) ^ 0:^^(0) [321. ^ 
is thus noteworthy that if one employs a more sophis- 
ticated, nonperturbative DSE truncation (37l. [38|. some 
of the infrared strength in the gap equation's kernel is 
shifted from Q into the dressed-quark-gluon vertex. This 
cannot materially affect the net infrared strength re- 
quired to explain obscrvablcs but does reduce the amount 
attributed to the effective coupling. 

We also used Eq. (|5|) to compute the masses of selected 
J = 0, 1 radial excitations and exotics, with the results 
presented in Table [TTJ The last column in the Table was 
prepared as follows. Wc fit the entries in each row to 
both m{uj) — constant and 

m (w) = uj(cq + Cioj), (11) 

then compute the standard-deviation of the relative error 
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0.4 


0.5 


0.6 


(J20 




0.214 




0.155 


0.147 


0.83 


rriQ— 


0.814 




0.940 


1.053 


0.03 


m 7T1 


1.119 




1.283 


1.411 


0.02 


m a 


0.970 




0.923 


0.913 


1.25 


m +- 


1.186 




1.252 


1.323 


0.34 


m lyi 


1.358 




1.489 


1.575 


0.14 


m p 


1.088 




1.046 


1.029 


1.22 


TOi-+ 


1.234 




1.277 


1.318 


0.60 


m pi 


1.253 




1.260 


1.303 


0.03 



TABLE II. Masses obtained with Eq. (jSJ) and Duj = 
(1.1 GeV) 3 . The subscript "1" indicates first radial excita- 
tion. The last column measures sensitivity to variations in 
r c := 020 «C 1 indicates strong sensitivity; and 020 ~ 1, 

immaterial sensitivity. Dimensioned quantities are reported 
in GeV. 

in each fit, <To for the constant and 02 for Eq. (|11[) . and fi- 
nally form the ratio: 020 = o^/co- In preparing the table 
we used Duj = (1.1 GeV) 3 . This has the effect of inflat- 
ing the 7r- and p-meson ground-state masses to a point 
wherefrom corrections to rainbow-ladder truncation can 
plausibly return them to the observed values [22j| . In 
this connection it is notable that the value reported for 
m a in Table [TT] matches estimates for the mass of the 
dressed-quark-core component of the cr-meson obtained 
using unitarised chiral perturbation theory [45j |. 

We have seen that with Eq. (JSJ) ground-state masses 
of light-quark pseudoscalar- and vector-mesons are quite 
insensitive to w e [0.4, 0.6] GeV. Any minor variation is 
described by a decreasing function. We emphasise that, 
for reasons which are understood [37], [H| , the behaviour 
of scalar- and axial-vector mesons is not well described 
in the rainbow-ladder truncation and hence the results 
reported herein do not allow reliable conclusions to be 
drawn about the w-dependence of their masses. 

In the case of exotics and radial excitations, the vari- 
ation with uj is described by an increasing function and 
the variation is usually significant. (That is also the case 
with Eq. © [USSl].) This is rcadil y understood. The 
quantity r u := l/ui is a length-scale that measures the 
range over which C/ir acts. For uj = 0; i.e., the i5 4 -function 
case, this range is infinite, but it decreases with increas- 
ing uj. One expects exotic- and excited-states to be more 
sensitive to long-range features of the interaction than 
ground-states and, additionally, that their masses should 
increase if the magnitude and range of the strong piece 
of the interaction is reduced because there is less binding- 
energy. (Recall cvflx(O) is a decreasing function of uj.) 

TablcHIlconfirms a known flaw with the rainbow-ladder 
truncation; viz., whilst it binds in exotic channels, it 
produces masses that are too light, just as it does for 
axial-vector mesons. Plainly, no predictions for exotics 
can be considered reliable unless the same formulation 
produces realistic predictions for axial-vector masses. It 



is similarly noticeable that m Xl is far more sensitive to 
variations in uj than is m pi ; and although m 1Tl < m Pl for 
uj = 0.4 GeV, the ordering is rapidly reversed. Thus, in 
conflict with experiment, one usually finds m 7ri > m pi in 
rainbow-ladder truncation. This, too, is a property of the 
truncation, which is insensitive to the details of G(k 2 ); 
e.g., the same ordering is obtained with a momentum- 
independent interaction [48| . 

It is probable these failings can be explained by the ac- 
tion of material, essentially nonperturbative corrections 
to the rainbow-ladder truncation which amplify effects 
that in quantum-mechanics would be described as spin- 
orbit interactions |37|. For example, the p„-mesons pos- 
sess nonzero magnetic- and quadrupolc-moments. This 
suggests that there is appreciably more dressed-quark or- 
bital angular momentum within these states than within 
7r n -mesons. Spin-orbit repulsion could significantly boost 
m pi and thereby produce the correct level ordering. In 
addition, the quark model describes axial-vector mesons 
as P-wave states and it has already been established that 
the non per turbative kernel corrections greatly boost their 
masses [38[ . Finally, exotic states appear as poles in ver- 
tices generated by interpolating fields with odd "time- 
parity" [49[ , a feature which magnifies the importance of 
orbital angular momentum within these states. 

We explored the efficacy of a form for the rainbow- 
ladder kernel at spacelike momenta which is a bounded, 
regular function that achieves its maximum value at 
k 2 = 0; viz., Eqs. (|3|), ([8]). In every respect tested, this 
form produces results for hadron observables that are 
at least equal to the best otherwise obtained in a com- 
parable approach. Given that the form we proposed is 
consonant with contemporary DSE- and lattice-QCD re- 
sults on the nature of the gap equation's kernel, in fu- 
ture studies it should replace other Ansdtze that fail in 
this respect. Justified, too, is repetition of numerous 
extant calculations; in particular, those relating to fea- 
tures of the phase transition at nonzero temperature and 
chemical potential, such as the novel effects discussed 
in Refs. [50, [H|, and properties of excited states other 
than their masses. In the latter connection, there is cer- 
tainly much to add. Herein, however, we will only report 
that, with Duj = (1.1 GeV) 3 and our favoured value of 
w = 0.6 GeV, 

/ Pl = -0.45/ p . (12) 

The sign is correct and significant, MS Ml* gued in Ref . 1461 . 

More important than repeating calculations, however, 
is movement beyond the rainbow-ladder truncation and 
its incremental improvement. Our results emphasise 
again that, within this circumscribed horizon, the out- 
standing questions in hadron spectroscopy, for example, 
do not have an answer. In order to use extant and forth- 
coming data as a tool with which to constrain the nature 
of QCD's /3-function, it is necessary to employ kernels 
that can more realistically handle systems other than 
ground-state vector- and flavoured-pscudoscalar-mcsons. 
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